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We consider Fermions in two dimensions with an attractive 
interaction in the singlet d-wave channel of arbitrary strength. 
By means of a Hubbard-Stratonovich transformation a statis- 
tical Ginzburg-Landau theory is derived, which describes the 
smooth crossover from a weak-coupling BCS superconductor 
to a condensate of composite Bosons. Adjusting the inter- 
action strength to the observed slope of H C2 at T c in the 
optimally doped high-T c compounds YBCO and BSCCO, we 
determine the associated values of the Ginzburg-Landau cor- 
relation length £ and the London penetration depth A. The 
resulting dimensionless ratio fcF^(O) ~ 5—8 and the Ginzburg- 
Landau parameter k = A/£ ~ 90 — 100 agree well with the 
experimentally observed values. These parameters indicate 
that the optimally doped materials are still on the weak cou- 
pling side of the crossover to a Bose regime. 
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I. INTRODUCTION 

The problem of the crossover from a BCS-like supec-. 
conducting state to a Bose condensate of local pairs&rEl 
has gained new interest in the context of high-Tc super- 
conductors. While there is still no quantitative micro- 
scopic theory of how superconductivity arises from dop- 
ing the, antiferromagnetic and insulating parent com- 
poundafl, it is clear that the superconducting state can be 
described in terms of a generalized pairing picture. The 
many body ground state is thus a coherent superposition 
of two particle states built from spin singlets in a rela- 
tive d-wave configurations. The short coherence length 
£ (0) ~ 10 — 20 A parallel to the basic Cu02-planes which 
is of the same order than the average interparticle spac- 
ing kp , however indicates that neither the BCS picture 
of highly overlapping pairs nor a description in terms of 
composite Bosons is applicable here. It is therefore of 
considerable interest to develop a theory, which is able 
to cover the whole regime between weak and strong cou- 
pling in a unified manner. On a phenomcnological level 
such a description is provided by the Ginzburg-Landau 
(GL) theory. Indeed, it is a nonvanishing expectation 
value of the complex order parameter ip which signals 
the breaking of gauge invariance as the basic characteris- 
tic of the superconducting state, irrespective of whether 
the pair size is much larger or smaller than the inter- 
particle spacing. A GL description of the BCS to Bose 
crossover was developed for the s-wave case by Drech- 
sler and one of the present authorscl in two dimensions 
and by Sa de Melo, Randeria and EngelbrechtEl for the 



three-dimensional case. In our present work the theory 
for the two-dimensional case is reconsidered, including 
a discussion of the Nelson-Kosterlitz jump of the order 
parameter and the generalization to the experimentally 
relevant situation of d-wave superconductivity. Moreover 
we also calculate the characteristic lengths £ and A and 
compare our results with measured properties of high-T c 
compounds. 

The remarkable success with which the standard BCS 
model has been applied to conventional superconductors 
relies on the fact that in the weak coupling limit the de- 
tails of the attractive interaction are irrelevant. By an 
appropriate rescaling of the parameters the properties of 
all weak coupling superconductors are therefore univer- 
sal. It is one of our aims here to investigate to which 
extent such a simplifying description also exists in more 
strongly coupled superconductors. Starting from a mi- 
croscopic model with an instantaneous attractive inter- 
action, we find that the resulting GL functional takes the 
standard form for arbitrary strength of the coupling. By 
adjusting a single dimensionless parameter to the mea- 
sured upper critical field near T c , we obtain consistent 
values for both the dimensionless ratio fcp£(0) » 5 — 8 
between the coherence length and interparticle spacing as 
well as the observed value n w 90— 100 of the GL param- 
eter in optimally doped high-T c compounds. Therefore, 
in spite of the rather crude nature of the original micro- 
scopic model, our GL theory is quantitatively applicable 
to strongly coupled superconductors which are far from 
the standard weak coupling limit, although not yet in the 
crossover regime to Bosc-like behaviour. 

The plan of the paper is as follows: in section II we 
introduce our microscopic model which has an attrac- 
tive interaction in the singlet d-wave channel. From this 
we derive, via a Hubbard-Stratonovich transformation, a 
statistical GL theory which is valid near T c . The relevant 
coefficients of the GL functional are calculated for arbi- 
trary interaction strength. In section III we discuss the 
appropriate microscopic definition of the order parame- 
ter and the evolution from the BCS to the Bose limit of 
the well known Kosterlitz-Thouless jump in the super- 
fluid density of two-dimensional superconductors. Ignor- 
ing the subtleties of the Kosterlitz-Thouless transition, in 
section IV we use a Gaussian approximation to determine 
the critical temperature and the associated value of the 
chemical potential at the transition. Finally in section V 
we determine the characteristic lengths £ and A near the 
transition for arbitrary strength of the coupling. Adjust- 
ing the coupling to the experimental values of the slope 
of H c2 near T c , we then determine the associated dimen- 
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sionless ratios kp£(0) and k — £/A. They agree rather 
well with the observed values in YBCO and BSCCO. 
A brief conclusion and a discussion of open problems is 
given in section VI. 



to a <-> —a.). For the derivation of a GL functional 
below, it is convenient to introduce pair operators b q via 



V k+q/2 C -k +1 C k+q -1 



(6) 



II. MICROSCOPIC DERIVATION OF THE GL 
FUNCTIONAL 

As a general model describing Fermions with an in- 
staneous pairwise interaction V(r) in a translationally 
invariant system, we start from the Hamiltonian (V is 
the volume of the system.) 
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with an arbitrary single particle energy e k which we 
will later replace by an effective mass approximation 
e k = ti 2 k 2 /2m. In the two-dimensional case, which we 
consider throughout, the Fourier transform V k _ k> of the 
interaction potential may be expanded in its relative an- 
gular momentum contributions I by 



V, 



fc-fc' 



V (k,k') 



1=1 



cos (7 ip)V l {k,k') 



(2) 



with <p the angle between k and k'. In the following 
we are only interested in d-wave pairs with symmetry 
d x 2_ y 2. We therefore omit all contributions I ^ 2 and 
also neglect the dependence on the absolute values k and 
k! of the momenta. Assuming the interaction is separa- 
ble, we thus approximate 



k-k' 



with 



-9 V k+q/2 V k'+q/2 
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and g a negative constant characterizing the strength 
of the attractive interaction. Finally the restriction to 
singlet pairing is incorporated trivially by considering 
only interactions between Fermions with opposite spins 
a' = —a. In this manner we obtain a Gorkov-like reduced 
interaction Hamiltonian 



H ' ~ 2V V k+q/2 V k'+q/2 



kk' q " 

C fc+q a C -fc -a C -k' -a C k'+q<j (^) 

(Note that the shift by q/2 in Eq. (|J) is necesssary to 
guarantee that the interaction is symmetric with respect 



The contribution H' may then be written in the form 

H ' = -T?T. b \K (?) 



of an attractive interaction between pairs of Fermions 
with opposite spin and total momentum q. In the follow- 
ing we want to derive a functional integral representation 
of the grand partition function 



Z = Tr e-^ H ~^ 



(8) 



which gives the standard GL theory as its mean field 
limit. Since we are interested in a superconducting state 
with a nonzero anomalous average (b q ) ^ 0, it is conve- 
nient to formally linearize the interaction term H' by a 
Hubbard-Stratonovich transformatiorfl. The grand par- 
tition function is thus expressed in terms of a functional 
integral 

Z = Jd 2 z exp(- ^- J dT5>(<7,T)| 2 ) L[z] (9) 

q 

over a complex valued c- number field z(q, r). Here 

0H 

L[z] = Tr Texp ( - i J dr H z (t)^ (10) 


is a functional of the auxiliary field z(q,r), which acts 
as a space- and 'time'-dependent external potential on a 
noninteracting Fermi system with Hamiltonian 

H Z {t) = ~ V) C ka C ka 

±J2(z(l,T)bl + z*(q,T)b q ). (11) 
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The physical interpretation of the c-number field z(g, r) 
is obtained by noting that its expectation value 



(z(q,T)) =-9(b q (r)) 



(12) 



is directly proportional to the anomalous average (6 q (r)). 
Thus up to some normalization constant, which will be 
determined below, the field z(q, r) is just the spatial 
Fourier transform of the complex order parameter ip(r,r) 
describing the superconducting state. It depends both on 
position and imaginary time r G [0, f3h] which is charac- 
teristic for a quantum GL functional. Since the d x 2_ y 2- 
symmetry in a rotational invariant system is comjtected 
with a one dimensional irreducible representationB, the 
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order parameter is still a simple complex scalar, similar 
to the more familiar isotropic s-wave case. 

Obviously the trace in the time ordered exponential 
in Eq. (10) cannot be calculated exactly. However it is 
straightforward to evaluate L[z] perturbatively in z. The 
naive justification for this is that close to T c the order 
parameter is small. Strictly speaking however, the func- 
tional integral in Eq. (||) requires to integrate over arbi- 
trary realizations of z(q, r). In order to obtain the stan- 
dard form of the statistical GL functional, however, the 
expansion is truncated at fourth order in the exponent of 
L[z]. In the language of field theory we are therefore cal- 
culating the bare coupling constants, which serve as the 
starting point for treating the behaviour at long wave* 
lcnghts. By a straightforward perturbative calculationE2l 
up to fourth order in z, the functional L[z] turns out to 
be of the form 

L[z] = Z Q exp ^X!E a ( 9 ' w "' \z(q,u}n)\ 2 



and replace 



l — £ 6(L 2 > 3) z(l)z*(2)*(3)**(l -2 + 3) 
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(13) 



Here Zq is the grand partition function of noninteracting 
Fermions while 



an 

z(q,uj n ) = J ' dr e w - T z{q,T) 




(14) 



is the Fourier transform of the r-dependence of z(q, r) 
with bosonic Matsubara frequencies uj n — 2im/ ((3K), n 
integer. In the quartic term we have used the short hand 
notation 1 = (qx,u>i), etc. The functions a and 6 can 
be expressed in terms of the normal state Green function 
(£ k = e k -H, u> n = 2w(2n + l)/(f3h)) 



G (k,uj ri 



(15) 



via 



h £j n 

■ Go(k + q i ui n + u} n )G (-k,-uj n ) (16) 

and a similar expression with four factors Go for 6. In 
order to obtain the standard form of a quantum GL func- 
tional, the coefficients a(q,uj n ) and 6(1,2,3) have to be 
expanded for small q and u) n . To lowest order in the spa- 
tial and temporal gradients of the order parameter, it is 
sufficient to keep only the leading terms in 



1 . h 2 q 2 

a(q, u) n ) — a(q, w n ) = a + c id fku ri 

g 2m 



(17) 



6(1,2,3) -> 6(0,0,0) 



(18) 



by its constant value at zero momentum and frequency. 
This expansion is valid, provided the contributions of or- 
der q 4 and ui 2 in Eq. (1^) are negligible. From an ex- 
plicit calculation of these higher order terms it may be 
shownlliJ that the order parameter must vary slowly on 
length scales of order ^ w K/ \fmEh with Eh the two 
particle binding energy introduced in Eq. ( p3| ) below. 
Physically, the length is just the radius of a bound 
state with energy In the weak coupling limit this 
length coincides with the standard BCS coherence length 
£o ~ frvF/T c which, for the clean limit considered here, 
is identical with the GL coherence length £(0) as defined 
in (^0|) . The standard form of the GL functional with a 
gradient term |V')/'| 2 is therefore valid provided the order 
parameter varies on scales larger than « £(0). With 
increasing strength of the coupling Eb, the pair radius 
£b decreases and thus the validity of the expansion (17) 
extends to variations on shorter length scales. Regard- 
ing the dependence on r, the requirement is that ip(r,r) 
must vary slowly on time scales r b sa h/E b . For weak 
coupling this is a rather large scale of order htpjT 2 (we 
set k,B = 1). Similar to the spatial dependence, how- 
ever, the necessary scale for the r-dependence of the or- 
der parameter for which the leading terms kept in (17) 
are sufficient, decreases with increasing coupling. Thus, 
in the Bose limit, the description of the time dependence 
of the order parameter by a first order deijiMative like 
in the well known Gross-Pitaevskii equationO becomes 
exact (see also section VI). 

With these approximations our GL functional in 



Z = Zq J D 2 z exp(- 



(19) 



finally takes the form 



v o 

+ dhz*(r,T)d T z(r, T ) + h - \z(r,r)A (20) 

which reduces to the familiar expression if z is indepen- 
dent of r. The coefficients a and 6 are given by 

a ^l-^^^ tanh f^ (21) 
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and 



J_ ^fc 



(22) 



Now the sum over wavevectors k in Eq. ( pi| ) diverges 
at large k and thus the bare value of a is undefined. In 
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the weak coupling limit this divergence is usually elimi- 
nated by argueing that the interaction is finite only in a 
thin shell around the Fermi surface. In the present case 
however, where the condensation in the strong coupling 
limit really affects the whole Fermi sphere, such a proce- 
dure is no longer passible. Instead, as was pointed out 
by Randeria et al.ll 2 ], we have to connect the bare cou- 
pling constant to the low energy limit of the two-body 
scattering problem. In two dimensions this relation is of 
the form 



and 



f 



4ttS 2 V E b J ' 
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(23) 



where £a - * oo is a high energy cutoff which precisely 
cancels the large k divergence on the right hand side of 
Eq. (pl|). The parameter Eb > is the binding energy 
of the two particle bound state in vacuum, which in fact 
must be finite in order to obtain a superconducting in- 
stability in two dimensions. In our present model, which 
neglects the dependence of V k _ k , on the absolute values 
of k and k', the existence of a bound state is indeed a nec- 
essary condition for superconductivity even in the case of 
d-wave pairing,palthough quite generally it only applies in 
the s-wave casetj. Since in the effective mass approxima- 
tion ej; = h 2 k 2 /2m which we are using throughout, the 
free Fermion density of states in two dimensions is con- 
stant, the coefficient a can now be calculated analytically 
in terms of Eb as 



a = y H ~ ] 6{jl) + H — ] + ln{ — ] 

^ gggfr) t ln^/2) \ 
1 2 ' 2 J ^ cosh 2 (/3£/2) ) V ' 



(7 = 0.577. . . is the Euler constant.). The coefficient b 
in (22) is finite without a cutoff and given by 



16nh 2 \ 2vr 2 



+ sgn(^ c ) J d£ 



tanh(/U/2) 



(25) 



Here 0(x) and sgn(x) are the well known Heaviside and 
sign functions. Moreover we we have replaced /3 and /1 
by their values at the critical point. Similarly, the values 
of the two remaining coefficients c and d at the critical 
point are 



8irh 2 I 
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00 
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tanh(/3 c ^/2) 
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(26) 



d = 



to tanh(/3 c /x c /2) 



8i:h 2 



(27) 



All four GL coefficients can thus be expressed essentially 
in analytical form for arbitrary strength of the interac- 
tion n Comparing these results with those for the s-wave 
caseQ, it turns out that up to a geometrical factor 3/2 in 
b the coefficients are identical at given values of j3 c and 
/i c , provided the two-particle binding energy Eb is simply 
identified with the corresponding s-wave value. 



III. ORDER PARAMETER AND 
NELSON-KOSTERLITZ JUMP 

In order to relate the formal auxiliary field z(r, r) in 
the functional ( p0| ) to the usual superconducting order 
parameter ip(r,T), the standard procedure in weak cou- 
pling is to take ^bcs = v2cz, which gives the conven- 
tional coefficient fi 2 /4m in front of |V?/> BCS | 2 . The gra- 
dient term is thus identical with the kinetic energy of a 
Schrodinger field for a single quantum mechanical parti- 
cle with mass to* = 2m, describing a pair built from coa- 
stituents with mass to. As pointed out by de GennesEil, 
however, the value of to* is arbitrary in principle, as long 
as one is considering the classical GL functional with ip 
independent of r. Indeed all measurable quantities ob- 
tained from the classical GL functional depend only on 
ratios like \ip\ 2 /m*. An arbitrary choice for to* can there- 
fore always be compensated by an appropriate rescaling 
of ip. This situation is changed, however, in a quan- 
tum mechanical treatment, where the order parameter 
also depends on r, i.e. dynamics enters. In this case 
there is a different natural normalization of ip = z in 
the coefficient of the ip* 9 T -0-contribution is just 
With this choice of normalization, the order pa- 



rameter ip(r,r) is precisely the c-number field in a coher- 
ent state path integralll3 associated with a genuine Bose 
field operator ip(r) with canonical commutation relations 
[ip(r), ip'(r')] = 5(r — r'). While this normalization is ev- 
idently the most appropriate one in the strong coupling 
Bose limit - where it agrees with the standard choice as 
we will see - it can be used for arbitrary coupling, even 
in the BCS- limit. Including the charge — 2e of a pair 
by generalizing the gradient to a covariant derivative in 
the standard way and adding the energy associated with 
the magnetic field h = V x A, the resulting free energy 
functional reads (cq is the velocity of light) 



F[r/>,A] 



, \- fd 2 r f dT (-^\iP\ 2 + -L- f^V + ^Ah 

ih J J V 2m* V i cq J 

V 

+ + f H 4 )+^- (28) 
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With this normalization, the three remaining indepen- 
dent coefficients now have a very direct physical inter- 
pretationli-3: the coefficient //* = — a/d is the effective 
chemical potential of the Bosons, m* = m d/c their ef- 
fective mass and g± — b/d 2 a measure of the repulsive 
interaction between the composite Bosons. 

In the following we will concentrate on the effective 
mass m* (at T c ) which, according to ( |26| , |27| ) is completely 
determined by the ratio T c /fx c . Now in the weak coupling 
limit \i c is equal to the Fermi energy (see section IV be- 
low) and thus m^/m vanishes like (T c /ep) 2 in the case 
of a BCS-superconductor. By contrast, for strong cou- 
pling fi c approaches ~Eh/2, i.e. is large and negative. In 
the Bose limit m* is therefore equal to 2m as expected. 
Relating T c //i c to the dimensionless coupling strength 
In E\y/ ep by the Gaussian approximation discussed in the 
following section, we obtain a monotonic increase of m* 
from exponentially small values to 2m as a function of 
coupling, as is shown in Fig. [|. 
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FIG. 1. The effective mass m+fm of the composite Bosons 
versus the binding energy Eb/tF- 

As was pointed out above, the mass m* of a Cooper 
pair defined in such a way cannot be observed in any 
static measurement like the penetration depth. To dis- 
cuss this, we consider the two dimensional current density 



-cq- 



5A(r) 



-2e|Vr — V<5 
m* 



which follows from (|28|) for a r-independent order pa- 
rameter ip = l^le^. For a spatially constant magnitude 
-0oo of the order parameter, this leads immediately to the 
London equation 



V x j(r) 



4e 2 |^ 



m*co 



-h 



(30) 



As was noted above, it is only the ratio \ip\ 2 /m+ which 
enters here and thus static magnetic properties are inde- 
pendent of the choice for m 4 . Specifically we consider a 
thin superconducting film with thickness S. The in-plane 



penetration depth A is then related to \tp c 
2_A2 



ia@ 



L„ = 



m*CQ 



(31) 



Here we have introduced a further length L s which is the 
effective magnetic penetration depth in a thin film. Typ- 
ically this length is of the order of one centimeter and 
thus for sample sizes which are smaller than that, mag- 
netic screening may be neglected. In such a situation the 
difference between a charged and a neutral superfluid be- 
comes irrelevant. A supercond-apting film thus exhibits a 
Kosterlitz-Thouless transitionO, in which the renormal- 
ized helicity modulus j(T) — h \ip\ (T)/m* jumps from 
2T c /n to zero at T c . Using (31) this jump translates inta 
one for the two dimensional screening length L s of sizeE£l 



(L S {T) T)\ 



1.96 K cm 



(32) 



where (f>o — hca/2e is the standard flux quantum. Con- 
sistent with our remarks above, this jump is completely 
universal and independent of m*, applying both to BCS- 
or Bose-like superconductors, provided L s is largjet; than 
the sample size and the density of vortices is lowtH. 

In order to define a proper superfluid density n s , we 
consider the relation between the order parameter ib and 
the microscopic anomalous average. From Eq. ( |l2|) we 
have ^bcs = —gV^c (b) ■ Neglecting the internal d-wave 
structure of the order parameter and the logarithmic fac- 
tor in (Eq) , it is straightforward to see that 



^Bcs(r-) w &(^+i(r) ip-i(r)) 



(33) 



with ipcr(r), a = ±1 the Fermionic field operators (The 
factor \n{2tk/ Eh) which is omitted in (|33| ) diverges as 
£a — * oo. This is a reflection of the fact that the product 
of two field operators at the same point can properly de- 
fined only with a cutoff.). Since is the radius of a pair, 
the relation j3^ ) shows that |V'bcs( t *)| 2 is just the areal 
density of pairs. This remains true even in the Bose limit 
where ^ while the product \jj + \{r) ^_i(r) eventually 
behaves like a true Bose field operator ip(r). The stan- 
dard definition n s = 2|-!/>bcs| 2 of the superfluid density 
can thus be applied for arbitrary coupling. By contrast, 
the Bose order parameter tji = y^m^/2m-0Bcs coincides 
with V'bcs only in the strong coupling limit. For weak 
coupling it is given by an expression like ( |33| ) but with the 
interparticle spacing kp 1 instead of as the prefactor. 
Thus ^(r)! 2 is essentially the probability density for two 
Fermions with opposite spin at the same point. In the 
BCS limit this density is exponentially small due to the 
large size of a pair. The superfluid density in turn is still 
of order one even in weak coupling and indeed at zero 
temperature n s must be equal to the full density n for 
any superfluid ground state in a translational system as 
considered herai3. Using n s = 2|^ B cs( r )| 2 j the relation 
( |32| ) can be rewritten in terms of a jump 

n ^ T c) = ^2 T c (34) 
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of the renormalized superfluid density. The superfluid 
fraction n s /n therefore has a jump of order T c /T^° se . 
Since this ratio approaches zero in weak coupling, there is 
a smooth crossover between the universal jump of the su- 
perfluid density in a Bose superfluidlij and the behaviour 
in a strict BCS model where n s /n = 2(1 — T/T c ) vanishes 
continuously near T c even in two dimensions. Indeed the 
BCS-Hamiltonian is equivalent to a model with an infi- 
nite range interaction, of strength V -1 for which mean 
field theory is exactE]. In the following we will neglect 
the subtleties associated with the Kosterlitz-Thouless na- 
ture of the transition, which is anyway masked by the 
coupling between different Cu02-planes in real high-T c 
superconductora^giving a three dimensional critical be- 
haviour near T^3. 



IV. GAUSSIAN APPROXIMATION 

In order to calculate directly observable quantities 
from our GL functional (p8|), we have to determine both 
the critical temperature T c and the corresponding chemi- 
cal potential fi c in terms of the binding energy Now 
it is obvious that an exact evaluation of the functional 
integral over ip(r, r) is impossible. We will therefore use 
the Gaussian approximation above T c , which is obtained 
by simply omitting the |?/>| 4 -term. With this approxima- 
tion our complete grand canonical potential f2 per volume 
takes the form 



diverges^. It is a straightforward generalization of the 



2m* 



- i hw n 



(35) 



The critical temperature and chemical potential then fol- 
low from the standard condition 



/x*(T c ,^ c ) = 



(36) 



for a bifurcation to a nonzero order parameter, and the 
particle number equation 



on 

n = -— = n a 

OjJL 



d 2 q 1 



Here 



n 



m ln(l + exp(/3/x)) 



irh 2 



(37) 



(38) 



is the number density of a free Fermion gas in two dimetu 
sions. Eq. (|3^) is identical with the Thouless criterionEa 
for a superconducting instability, which is equivalent to 
the condition that the ladder approximation to the exact 
pair field susceptibility 



Xpair = J dr (6 (r) & o (0)) 
o 



(39) 



usual gap ec 
equation 
have 



uation to arbitrary coupling. The number 
deserves some more comments. Since we 



(40) 



quite generally, it is easy to see that <9 M /i* = 2 at — (3 C 
and // = fi c . Therefore Eq. (^) has the simple intu- 
itive interpretation that the total number of particles is 
split into the number of free Fermions still present at 
(/3 c ,/z c ) plus the number of Fermions already bound to- 
gether in pairs, whose mean occupation number is just 
the Bose distribution. Now formally this distribution 
function arises from the summation over the Matsubara 
frequencies uj n in Eq. (|3^) precisely because our coeffi- 
cient a(q,Lu n ) has been expanded only to linear order in 
cu n . The omission of the higher order terms in this ex- 
pansion is therefore connected with neglecting scattering 
state contributions, which would give an additional term 
in Eq. (^) beyond the completely free and fully bound 
number of Fermions. Such a contribution is important 
in the three-dimensional case where true bound st 
exist only beyond a critical strength of the coupling 
For our present discussion of the problem in two di- 
mensions, however, there are only free Fermions or true 
bound states. Therefore there is no contribution in Eq. 
(|37| ) from scattering states and one expects that the ex- 
pansion of a(q,cj„) to linear order in uj n is reliable at 
arbitrary strength of interaction. 

There is however a rather different problem which ap- 
pears in the two-dimensional case. As was discussed 
above, a superconducting transition exists only in the 
Kosterlitz-Thouless sense. This problem shows up in our 
Gaussian approximation, since the Bose integral in Eq. 
(|37| ) diverges at [i+ = 0. Now at this level of approxima- 
tion this is just a reflection of the fact that T c — for an 
ideal Bose gas in two dimensions, because - as pointed 
out above - the omission of the l^^-term corresponds to 
neglecting the repulsive interaction between the Bosons. 
From our analytical results ( p5| ) and ( f27| ) for b and d, 
which do not contain Eb, it is however straightforward 
to calculate the effective interaction g± for arbitrary cou- 
pling. It turns out that g+ is a monotonically decreasing 
function of the coupling. In the limits (3 c fi c — » oo of a 
BCS-like or /3 c ^t c — > — oo of a Bose-like system, we find 



BCS 
Bose 



(41) 



Thus, in two dimensions, there is always a finite repul- 
sive interaction-between the pairs, which is of purely sta- 
tistical origirfit-j. In particular g+ remains finite in the 
Bose limit, where it arises from processes with a virtual 
exchange of one of thej-cpnstituent Fermions in a Bose- 
Bose scattering processed. The fact that is very large 
in the weak coupling limit is simply a consequence of 



G 



the large pair size £o ~ Pc^f (note that /i c = ef in the 
weak coupling limit), but does not imply that the |V'| 4 - 
term is particularly relevant in this regime. On the con- 
trary, using the .-Gaussian approximation, it is straight- 
forward to showed that in this limit the product g*(|i/>| 2 ) 
which effectively renormalizes the Boson chemical poten- 
tial /Lt* , is of order ^JE^tF which is roughly T c in the weak 
coupling limit. For BCS-like superconductors the 
contribution is therefore irrelevant except very close to 
T c , a fact which is well known from the standard theory 
of conventional superconductors. Now the finite value of 
g+ guarantees that even in two dimensions there is a finite 
critical temperature below which the superfluid density 
is nonvanishing. Unfortunately it is not possible to incor- 
porate the Kostcrlitz-Thouless nature of the transition in 
an approximate treatment of the GL functional. However 
considering the effectively three-dimensional structure of 
high-T c superconductors, this problem may be circum- 
vented by including the motion in the direction perpen- 
dicular to the planes. A very simple method to incorpo- 
rate this is provided by adding a transverse contribution 
ej_ to the kinetic energy bynJ 



h 2 q 2 h 2 q 2 h 2 q\ 
2m* 2m* 2mi 



(42) 



whose average is equal to the thermal energy (ej_) = T. 
Replacing the integral in Eq. (5^) by 



d 2 q 



L 



d 2 q f dq± 



(27T)* 



2tt 



(43) 



with L± — 2tt I ' yj (qf) we obtain an effectively three- 
dimensional system. This becomes evident by writing 
the contribution of the bound pairs in (fj^ ) in the form 



m* y/J3 
irk 2 



de 



exp[, 



+ e) 



(44) 



The density of states is thus proportional to y/e as in 
three dimensions, making n' finite at /i* = 0. Although 
rather crude, this approximation gives a value for T c 
in the Bose limit, which is very close to the Kjosterlitz- 
Thouless value for the transition temperatureO 



0.89 



nh 2 n 



(45) 



of a dilute hard core Bose gas on a latticed with Boson 
mass 2m and number density n/2. 

Using our results for_the GL coefficients a, c and d 



and the replacement (44) in the number equation, we 
can now determine both T c and /i c for arbitrary cou- 
pling from /it* = and Eq. (|3^). The corresponding 
results are shown in Figs. 2 and 3 in units of the char- 
acteristic energy scale and as a function of the di- 
mensionless effective coupling E^jtp. For weak coupling 
Eb <C £f, the critical temperature is monotonic in E^, 



behaving like T c y/JUb^F- For intermediate coupling, it 
exhibits a maximumQ. A similar but less pronounced be- 
haviour is found in .three dimensions, again in the Gaus- 
sian approidmationErQ. In a more refined self-consistent 
treatmentE3, however, the critical temperature is a mono- 
tonically increasing function of the coupling. It is likely 
that the same situation also applies in two dimensions, 
but unfortunately there is at present no quantitative the- 
ory taking into account the repulsive interaction between 
the pairs in this case. In the Bose limit the transition 
temperature becomes independent of the original attrac- 
tive interaction and is completely determined by the Bo- 
son density n/2 and the effective mass m* — > 2m. It 
approaches a value of about 10% of the Fermi energy, 
which is likely to be an upper limit for T c in the present 
problem. 




FIG. 2. The normalized critical temperature T c versus the 
dimensionless coupling strength. 

The chemical potential u c at the transition decreases 
monotonically from its weak coupling value cf to —Eb/2 
in the Bose limit. It changes sign at l^E^/cp) « —1.1, 
where the behaviour crosses over from BCS- to Bose-like 
(In Fig |^ we have supressed a small dip in fi c around 
these couplings which is an artefact of the pronounced 
maximum in T c ). Apart from the chemical potential /i c , 
the nature of the transition can also be inferred from 
evaluating the number of preformed Bosons at T c . This 
quantity, which is just half of the contribution ( f44| ) to 
the number equation is shown in Fig. ^. It is obvious 
that the nature of the phase transition changes rather 
quickly in a range of couplings between l^Eb/ep) ~ —2 
and —1. For smaller couplings, the density of preformed 
pairs near T c is essentially negligible and binding occurs 
simultaneously with condensation. On the other hand, 
for ln(i?b/ej?) > —1 basically all pairs are already present 
above T c and the transition to superconductivity is that 
of a true Bose system. 
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-4 -2 

Ln(E b /e F ) 

FIG. 3. The normalized chemical potential fi c at the criti- 
cal point versus coupling. 




-6 -4 -2 

Ln(E b /8 F ) 

FIG. 4. Contribution n' of bound pairs at the critical point 
to the total Fermion number density n. 

V. CHARACTERISTIC LENGTHS 

In the following we want to determine both the coher- 
ence length £ and the penetration depth A as a function 
of the coupling. The former is defined both above and 
below T c and may be obtained from the GL coefficients 
simply via 



2m a 



(46) 



The definition of a penetration depth in a two-di- 
mensional superconductor has been discussed in section 
III. Within the Gaussian approximation we may replace 
l^oo) 2 in (31) by n+/g+ which leads to 



A 2 = Ai 



bn 
icjoj 



(47) 



Here we have introduced the bare value Xl of the London 
penetration depth defined by 



A 2 = " lc o 



with ri3 = n/S the nominal three-dimensional carrier den- 
sity. Note that both £ and A have been written in terms 
of the original static GL-coefficients a, b and c, in order 
to stress that the characteristic lengths are independent 
of the normalization of ip, i.e. the kinetic coefficient d 
does not enter here. Since the coefficient a vanishes at 
the transition, both £ and A diverge. Now in a full treat- 
ment of the GL-functional, including the ip -term, the 
behaviour very close to T c in a single layer would be of the 
Kosterlitz-ThouLess type. The correlation length would 
thus diverge likeB ln£ ~ \T - T c |~ 1/2 while A" 2 would 
jump from zero to a finite value below T c (see (32)). In 
the three-dimensional case, the behaviour very close to 
T c is that of a 3d XY-mpdel with nontrivial but well 
known critical exponentsEa. In our Gaussian approxima- 
tion this complex structure is replaced by a simple mean 
field behaviour. However there is a subtle point even at 
this level of approximation. Indeed the coefficient a de- 
pends both on temperature and chemical potential and it 
is only in the strict BCS-limit, where the latter is a fixed 
constant equal to ep. With increasing coupling, however, 
the chemical potential changes and thus the relevant limit 
close to T c is to consider a(T,fi(T)) as T — > T c . Now 
by using the exact relation (TO), it is straightforward to 
showta that a(T,fj,(T)) ~ (T - T c ) 2 vanishes quadrati- 
cally near T c . The resulting critical exponent for the cor- 
relation length would thus be v = 1. Indeed this is the 
exponent expected for an ideal Bose gas in three dimen- 
sions, to which our Gaussian approximation is effectively 
equivalent. Now in order to allow a comparison of our 
results with measured values of £ and A, which are found 
to obey a mean field behaviour with v = 1/2 except very 
close to TjEB, we neglect the temperature dependence of 
the chemical potential near T c . As a result 



a(T, (j, c ) = a' 



T-T n 



(49) 



vanishes linearly near T c , giving the standard mean field 
divergence of £ and A. Obviously this approximation is 
only reliable on the weak coupling side of the transition. 
As we will see below, however, this is indeed the rele- 
vant regime even in high-T c superconductors. We thus 
expect that our results are at least qualitatively reliable 
for these systems. For strong coupling, the derivative 
a' of a(T,/j, c ) near T c vanishes like exp(— f3 c Ei,/2). As a 
result, the characteristic lengths £ and A would increase 
exponentially in the Bose limit which is unphysical how- 
ever. We have therefore restricted our calculation of the 
GL coherence length £(0) defined by 



2m a' 



(50) 



to coupling strengths lnffib/ep) smaller than about — 1, 
where the system starts to cross over to Bose like be- 
haviour. The corresponding result for £(0) in units of 



(48) kp 1 is shown in Fig. g. 
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Ln( E b /s F ) 

FIG. 5. The Ginzburg-Landau coherence length £(0) versus 
the coupling strength. 

It exhibits the expected decrease of the coherence 
length from its weak coupling limit 



k F m\ 



BCS 



7C(3) ££ 
8tt 2 T c 



(51) 



to values of order one near the crossover regime, before 
it starts to rise again. As was noted above our approxi- 
mations in determining £(0) are no longer reliable in this 
regime. While in three dimensions £(0) is expected to 
increase like 1/ \Jclb with &Brr= ^ a F —* the rele- 
vant Bose-Bose scattering lengthEl, the actual behaviour 
in two dimensions is unknown. Fortunately, however, this 
problem does not arise in determining the GL-parameter 
k = A/£ which, in two dimensions, can be expressed as 



/ m 
\8irh 2 



b5 
c 2 r c i 



1/2 



(52) 



Here we have introduced the equivalent of the classi- 
cal electron radius r c \ = e 2 /mcQ where m is the band 
mass. Since the problematic coefficient a has dropped 
out in k, we can use (|52] ) to determine the GL-parameter 
in the whole regime between weak and strong coupling. 
In the BCS-limit k is exponentially small, behaving like 
k bcs ~ T c /eF ■ {8/rd) 1 / 2 . The associated penetration 
depth A(0) = \/3/4 A^ is thus essentially equal to the 
London value A^. By contrast, for the Bose case k ap- 
proaches a constant value 



3<S\V2 



(53) 



which is large compared to one, since 6 3> r c \ for realistic 
values of the sheet thickness 5. The complete dependence 
of k in the crossover regime is shown in Fig. o. 



-6 -4 -2 

Ln( E b /e F ) 

FIG. 6. The GL parameter k versus coupling. 

It is a monotonic function of the binding energy. Thus 
with increasing coupling there is always a transition from 
type I to type II behaviour in two dimensions, even for 
a clean superconductor with no impurities as considered 
here. 

With these results, we are now in a position to com- 
pare our simple model with experimental values for op- 
timally doped high-T c superconductors. Since the criti- 
cal temperature is exponentially sensitive to the coupling 
strength and also is likely to be considerably reduced 
compared to our result in Fig. 2 by fluctuation effects in 
the crossover regime, we have refrained from taking T c 
as a reliable parameter for adjustment. Instead we have 
used the measured values of the slope of the upper criti- 
cal field H C2 near T c , which determines the GL coherence 
length £(0) via 



d\nT 



T=T C 



2< 2 (0) 



(54) 



In order to fix the dimcnsiotaJcss coupling strength Eb/eF 
from the measured valuesEJ = 14 A for optimally 
doped YBCO (T c = 92 K) &nM £(0) = 24 A for the cor- 
responding compound BSCCO (T c = 1Q7 K), we need 
the in plane carrier densities n which determine an effec- 
tive value of kp — \/2im (It should be pointed out that 
k F is introduced here only as a measure of the carrier 
density. In fact due to the strong attractive interaction, 
the actual momentum distribution of the Fermions just 
above T c will be far from the standard Fermi distribu- 
tion, except in the BCS limit.). The three-dimensional 
carxier densities in optimally doped YBCO and BSCCO 



ar< 



n 3 = 4 • ltr 1 cm" 



an 



dE3 Tla = 



rUJ cm re- 
spectively. With the corresponding valuesE3 § = 5.84 A 
and 5 = 9.27 A of the effective sheet thickness, the re- 
sulting Fermi momenta are kp — 0.38 A for YBCO 
and kp = 0.34 A for BSCCO. The dimensionless ra- 
tio fcp£(0) = 5 and 8 between the coherence length and 
the average interparticle spacing then allows us to deter- 
mine the effective coupling strength. From Fig. |we find 
that \n{Eb/e F ) is equal to —5.1 for YBCO and —6.0 for 
BSCCO. As Figs. || and ^ show, these coupling strengths 



describe superconductors which are still on the weak cou- 
pling side of the crossover from BCS to Bose behaviour. 
For instance the density of preformed Bosons at T c is 
less then 1% in both cases (see Fig. In order to check 
whether our description is consistent, we determine the 
associated values of the GL parameter n near T c . From 
the above values of 8 and the band massesa m — 3.5 m e 
for YBCO ancCJ m = 4.4 m e for BSCCO (m e is the free- 
electron mass.) we find that k Bosc is equal to 1.5- 10 3 and 
2.1 ■ 10 3 respectively for the two compounds considered 
here. From Fig. | we thus obtain k — 99 and k = 87 
for optimally doped YBCO and BSCCO. These numbers 
agree very well with, the experimcnt^tandetermined val- 



ues of k, which ar k cxd = 100 andf-3 



'•'cxp 



= 86. Our 

simple one parameter model therefore gives a consistent 
quantitative description of the characteristic lengths £ 
and A. 



VI. CONCLUSION 

To summarize, we have studied the crossover in the su- 
perconducting transition between BCS- and Bose-like be- 
haviour within a GL description. It has been found that 
optimally doped high-T c superconductors are still on the 
weak coupling side of this crossover, although they are 
certainly far away from the BCS-limit. Our microscopic 
model is characterized by a single dimcnsionlcss param- 
eter, similar to the familiar BCS-Hamiltonian. While 
the GL functional has the standard form for arbitrary 
coupling its coefficients depend strongly on the coupling 
strength. The crossover, at least in two dimensions, is 
essentially identical for the s- or d-wave case. For static 
properties, only two of the relevant coefficients a, b and c 
are independent, since the normalization of the order pa- 
rameter is arbitrary. Fixing £(0) from experiment there- 
fore leaves only one further parameter - for instance k - 
as an independent predicted quantity. The good agree- 
ment of k with measured values supports our conclusion 
that the optimally high-T c compounds are intermediate 
between BCS and Bose behaviour. Since the crossover 
regime is rather narrow, however, weak coupling theories 
are still a reasonable approximation for the relevant cou- 
pling strengths. This is consistent with the empirical fact 
that a weak coupling approach apparently works well in 
many cases. 

Evidently there are a number of important open ques- 
tions. They may be divided into two classes: the first 
one concerns the problem of a better and more complete 
treatment of our model itself. The second class is related 
to the problem, to which extent this model is applicable 
to high-T c compounds and what are the necessary in- 
gredients for a more realistic description. Regarding our 
microscopic Hamiltonian as a given model, it is obvious 
that our treatment of the associated GL phenomenology 
is still incomplete. In particular the Gaussian approxi- 
mation is obviously not sufficient to give a quantitatively 



reliable result for T c at intermediate coupling. Moreover, 
the behaviour of the characteristic lengths £ and A in the 
strong coupling limit is completely unknown. In order 
to go beyond the Gaussian approximation, it is neces- 
sary to include the pair interaction (i. e. the |i/'| 4 -term) 
properly. Progress in this direction hasJaeen made in 
the three-dimensional case by HaussmanncJ and very re- 
cently by Pistolesi and StrinatiEEl. Using a self-consistent 
and conserving approximation for the Green- and vertex 
functions, Haussmann obtained a smooth increase of T c 
with coupling, thus eliminating the unphysical maximum 
in the crossover regime found in the Gaussian approxima- 
tion. This approach is rather different from our present 
one and requires extensive numerical work. Pistolesi and 
Strinati have performed an essentially analytical calcu- 
lation of the correlation length £o at zero temperature, 
using a loop expansion in a functional approach simi- 
lar to our present one. They have shown that the pair 
radius £(, « h/\ZmE(, coincides with £ not only in the 
BCS-limit, but down to values around kp^Q « 10. Sim- 
ilar to our results for the GL coherence length £(0) in 
Fig. M, kp^Q reaches a minimum of order one in the 
crossover regime before it starts to rise again. Since 
in three dimensions the behaviour at strong coupling is 
that of a weakly interacting Bose gas with scattering 
length as — 2a-F ~ » OlJc-3, £o eventually increases like 
kp 1 1 ' \Jkpa,B while £j, rs clb approaches zeroQ. Unfortu- 
nately for the two-dimensional case, where the Bason in- 
teraction is finite even at very strong couplingErE-3, the 
Kosterlitz-Thouless nature of the transition makes it very 
difficult to improve upon the simple approximations used 
here. A first step in this direction was taken by Traventj. 
He showed that interactions between the pair fluctua- 
tions guarantee a nonvanishing superfluid density at fi- 
nite temperature, in agreement with our arguments be- 
low Eq. (41). However there seems to be no quantitative 
calculation of the coherence length in a two-dimensional 
Bose-like regime even near zero temperature. 

A different problem we want to mention here is that of 
the proper time dependent GL theory. Since our quan- 
tum GL functional is derived from a microscopic Hamil- 
tonian, in principle it contains the complete information 
about the order parameter dynamics, at least as far as 
intrinsic effects are concerned. Neglecting higher order 
terms in the expansion in u n , the resulting equation of 
motion for the order parameter in real time t isLO 



d ■ ih 



dip{r,t) _ 6F[il>] 



dt 



6ip*(r,t) 



(55) 



Due to the analytic continuation, the coefficient d — 
d% + «g?2 has now acquired a finite imaginary part c?2 > 
which describes irreversible relaxation. For a better com- 
parison with the standard literature, it is convenient to 
choose the conventional order parameter ipBGS = V2cz, 
where the prefactor of the gradient term is H 2 /4m. With 
this choice of normalization, the kinetic coefficient d\ — 
m* /2m is identical with the effective mass discussed in 
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section III. It is then evident from Fig. [j] that a Gross- 
Pitaevskii-like dynamics where d\ — 1, is only valid in 
the Bose limit, while d\ ~ (T c /ei?) 2 is exponentially small 
for weak coupling. Indeed for BCS-like systems it is c?2 
which is dominant, being of order T c jtp in the three- 
dimensional casecH. This result reflects the fact that for 
weak coupling superconductors the order parameter dy- 
namics is purely relaxing. Going beyond the BCS-limit, 
the associated kinetic coefficient di has only been evalu- 
ated in three dimensionsEO, where it exhibits a maximum 
at intermediate coupling. Its behaviour in the Bose-limit 
and in two dimensions in general, however, is completely 
unknown. Since the incorporation of scattering states in 
the three-dimensional case requires to go beyond the lin- 
ear expansion in to n , it is likely however, that a simple 
first order equation like Eq. (Bq) is in fact not appro- 
priate for describing the dynamics at intermediate cou- 
pling. It is only at very low temperatures where the sit- 
uation becomes simple again. Indeed from quite general 
argumentscJ, the dynamics as T — > is expected to be 
of the Gross-Pitaevskii form, irrespective of the strength 
of the coupling. Finally we mention that a proper micro- 
scopic calculation of the complex coefficient d is relevant 
for understanding the Hall effect in high-T c compoundsEa. 

Concerning the question to which extent our model 
is really applicable to high-T c superconductors, it is ob- 
vious that most of the complexity of these systems is 
neglected here. In particular we have assumed that the 
normal state is characterized by a given density of effec- 
tive mass Eermions with some instantaneous attractive 
interactionEj. Such a system will certainly be very dif- 
ferent from a conventional Fermi liquid for intermediate 
or strong coupling. Our conclusion that we are still on 
the weak coupling side of the crossover, with a negli- 
gible density of preformed pairs, is consistent with the 
phenomenology of optimally (and perhaps overdoped) 
high-X^ superconductors, however it is certainly inappro- 
priate for the underdoped cuprates. Indeed these sys- 
tems exhibit a gap far above T c which may be inter- 
preted in terms of preformed Bosons. A GL descrip- 
tion of underdoped compounds was vecy recently devel- 
oped by Geshkenbein, Ioffe and Larkinc9. Assuming that 
Bosons form far above T c only in parts of the Fermi sur- 
face and coexist with unpaired Fermions through T c , they 
obtain a reasonable description of the phenomenology of 
certain underdoped materials. This behaviour is quite 
different from that obtained in our model, which com- 
pletely neglects any effects of the Coulomb repulsion and 
band structure. It is obvious that for a quantitative de- 
scription of high-T c superconductors both Coulomb cor- 
relations and band structure effects have to be included, 
which requires to use lattice models. They provide a 
quantitative description of these complex systems at least 
in the normal state and allow the calculatiop-,of micro- 
scopic properties like spectral functions, etdiil. Unfor- 
tunately with these models it still seems impossible to 
really explain how d-wave superconductivity arises from 
the strongly spin and charge correlated normal state. As 



a result, effective models like the negative-U Hubbard 
model are often used-te, discuss microscopic properties of 
high-T c compoundsc3c3. Our present approach is more 
phenomenological, starting from a model in which all mi- 
croscopic details arc neglected except for the fact that we 
have a strong pairing interaction in a system of Fermions 
with given density. The advantage of such an approach is 
that it allows a simple calculation of the relevant lengths 
£ and A and quantities following from that like the criti- 
cal fields. The fact that the resulting GL theory gives a 
consistent description of optimally doped cuprates indi- 
cates that at least at this level the microscopic details are 
not relevant. Certainly our results supporting this view 
are quite limited so far and it is necessary to investigate 
this further. Since the coefficients of the GL functional 
near T c are quite generally determined by the properties 
in the normal state, an interesting future direction would 
be to see whether superconducting properties below T c 
can quantitatively be obtained from the GL functional 
by properly incorporating the anomalous behaviour in 
the normal state. 



ACKNOWLEDGMENTS 

One of the authors (W. Z.) would like to thank A. 
J. Lcggctt for his hospitality at the University of Illinois 
where this work was completed and for useful discussions. 
Part of this work was supported by a grant from the 
Deutsche Forschungsgemeinschaft (S. S.). 



1 A. J. Leggett, in Modern Trends in the Theory of Con- 
densed Matter, edited by A. Pekalski und J. Przystawa 
(Springer, Berlin, 1980). 

2 P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 
195 (1984). 

3 M. Randeria, in Bose-Einstein Condensation, edited by D. 
Snoke and S. Stringari (Cambridge University Press, 1994). 

4 D. J. Scalapino, Phys. Rep. 250, 329 (1995). 

5 J. Annett, N. Goldenfeld and A. J. Leggett, in Physical 
Properties of High Temperature Superconductors, Vol. 5, 
edited by D. M. Ginsberg (World Scientific, New Jersey, 
1996). 

6 M. Drechsler and W. Zwerger, Ann. der Phys. 1, 15 (1992). 

7 C. A. R. Sa de Melo, M. Randeria and J. R. Engelbrecht, 
Phys. Rev. Lett. 71, 3202 (1993). 

8 B. Miihlschlegel, in Path Integrals and their Applications in 
Quantum, Statistical, and Solid State Physics, edited by G. 
J. Papadopoulos und J. T. Devreese (Plenum Press, New 
York, 1977). 

9 M. Sigrist and K. Ueda, Rev. Mod. Phys. 63, 239 (1991). 
S. Stintzing, PhD thesis, Ludwig-Maximilians-Universitat 

Miinchen 1996, unpublished. 



11 



11 E. P. Gross, J. Math. Phys. 4, 195 (1963). 43 R. Micnas, M. H. Pedersen, S. Schafroth, T. Schneider, J. 

12 M. Randeria, J.-M. Duan and L.-Y. Shieh, Phys. Rev. B J. Rodriguez-Nunez and H. Beck, Phys. Rev. B 52, 16223 
41, 327 (1990). (1995). 

13 P. G. de Gennes, Superconductivity of Metals and Alloys 
(Addison- Wesley, Redwood City, 1989) 

14 W. Zwerger, in Proceeedings of the Fourth International 
Conference on Path Integrals from meV to MeV, edited 
by H. Grabert, A. Inomata, L. S. Schulman and U. Weiss 
(World Scientific, Singapore, 1992). 

L. S. Schulman, Techniques and Applications of Path Inte- 
gration (Wiley, New York, 1981), Chap. 27. 

16 B. I. Halperin and D. R. Nelson, J. Low Temp. Phys. 36, 
599 (1979). 

17 P. Minnhagen, Rev. Mod. Phys. 59, 1001 (1987). 

18 A. J. Leggett, Physica Fennica 8, 125 (1973). 

19 D. R. Nelson and J. M. Kosterlitz, Phys. Rev. Lett. 39, 
1201 (1977). 

20 B. Mublschlegel, J. Math. Phys. 3, 522 (1962). 

21 H. Keller and T. Schneider, Int. J. Mod. Phys. B 8, 487 
(1994). 

22 D. J. Thouless, Ann. of Phys. 10, 553 (1960). 

23 A. Tokumitu, K. Miyake and K. Yamada, Phys. Rev. B 47, 
11988 (1993). 

24 R. Haussmann, Z. Phys. B 91, 291 (1993). 

25 V. P. Gusynin, V. P. Loktev and I. A. Shovkovyi, Sov. 
Phys. JETP 80, 1111 (1995). 

26 In the continuum the result (^) is modified by a double 
logarithmic factor, which is neglected here however, see D. 
S. Fisher and P. C. Hohenberg, Phys. Rev. B 37, 4936 
(1987). 

27 R. Haussmann, Phys. Rev. B 49, 12975 (1994). 

28 S. Kamal, D. A. Bonn, N. Goldenfeld, P. J. Hirschfeld, R. 
Liang and W. N. Hardy, Phys. Rev. Lett. 73, 1845 (1994). 

29 X.-G. Wen and R. Kan, Phys. Rev. B 37, 595 (1988). 

30 U. Welp, W. K. Kwok, G. W. Crabtree, K. G. Vandervoort 
and J. Z. Liu, Phys. Rev. Lett. 62, 1908 (1989). 

31 I. Matsubara, H. Tanigawa, T. Ogura, H. Yamashita and 
M. Kinoshita, Phys. Rev. B 45, 7414 (1992). 

32 K. Semba, T. Ishii and A. Matsuda, Phys. Rev. Lett. 67, 
769 (1991). 

33 D. R. Harshman and A. P. Mills, Jr, Phys. Rev. B 45, 
10684 (1992). 

34 L. Krusin-Elbaum, A. P. Malozemoff, Y. Yeshurun, D. 
C. Cronemeyer and F. Holtzberg, Phys. Rev. B 39, 2936 
(1989). 

35 F. Pistolesi and G. C. Strinati, Phys. Rev. B 53, 15168 
(1996). 

36 S. V. Traven, Phys. Rev. Lett. 73, 3451 (1994). 

37 M. Stone, Int. J. Mod. Phys. B 9, 1359 (1995). 

38 For a phenomenological discussion see A. T. Dorsey, Phys. 
Rev. B 46, 8376 (1992). 

39 For a treatment of the strong coupling limit for systems 
with a retarded interaction in terms o f Eliashberg theory 
see G. Varelogiannis and L. Pietronero, :ond-mat/9507072| . 

40 V. B. Geshkenbein, L. B. Ioffe and A. I. Larkin, Phys. Rev. 
B 55, 3173 (1997). 

41 G. Dopf, J. Wagner, P. Dieterich, A. Muramatsu and W. 
Hanke, Phys. Rev. Lett. 68, 2082 (1992). 

42 M. Randeria, N. Trivedi, A. Moreo and R. T. Scalettar, 
Phys. Rev. Lett. 69, 2001 (1992). 



12 



